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Abstract. 

Recently, it was proposed that phonons are specularly reflected below ^0.5 K in ordinary 
single-crystal samples of high-T c cuprates, and that the low-temperature thermal conductivity 
k should be analyzed by fitting the data up to ~0.5 K with the formula k = aT + bT a where 
a is a fitting parameter. Such an analysis yields a result different from that obtained from the 
conventional analysis in which the fitting is restricted to a region below ~0.15 K, where k has 
been found to obey k = aT + bT 3 very well, which means that phonons are in the boundary 
scattering regime without specular reflection. Here we show that the proposed new analysis is 
most likely flawed, because the specular phonon reflection means that the phonon mean free 
path £ gets longer than the mean sample width, while the estimated £ is actually much shorter 
than the mean sample width above ~0.15 K. 



1. Introduction 

Because high-T c cuprate superconductors have a d-wave superconducting gap that has 
four "nodes" where the gap vanishes, there are low-energy quasiparticles (QPs) in the 
superconducting state of cuprates |lj. Those QPs, called nodal QPs, are extended and have 
well-defined wave vectors, and they can carry heat even at very low temperature. The thermal 
conductivity k of cuprates below T c has turned out to be a useful probe for studying the nature 
of the superconducting state, because the physics of the nodal QPs are directly reflected in k. In 
particular, at very low temperature where the QPs are in the residual elastic scattering regime, 
the electronic thermal conductivity K e is expected to be approximated by a "universal" thermal 
conductivity no, which is expressed as [21 [3] 

no _ /-up + V2_\ ^ fc| nvp_ ,^ 

where n is the number of Cu02 planes per unit cell, c is the c-axis lattice constant, and vp 
{y<i) is the QP velocity normal (tangential) to the Fermi surface at the gap node. This kq is 
called "universal" in the sense that it is expected to be independent of the QP scattering rate 
7; this happens because the effects of disorder to create QPs and to scatter them would exactly 
cancel each other in a certain range of 7 [T]. Thus, the bulk measurement of K e at low enough 
temperature would give us microscopic information about the gap. 



However, in metals both electrons and phonons carry heat, so k is a sum of K e and the phonon 
thermal conductivity n p . This causes a fundamental complication in the data analysis, that is, 
one has to separate K e from k p in the measured k. Such a separation can be accomplished 
if K e and k p have distinct and well-defined temperature dependences; actually, at low enough 
temperature, K e ~ aT is expected in the residual elastic scattering regime, and k p ~ 6T 3 is 
expected in the boundary scattering regime of phonons. Hence, experimentalists often employ 
dilution refrigerators to measure the thermal conductivity in the mK region, where the heat 
transport is usually found to be in the above-described regime. Indeed, in the case of cuprates, 
we have recently shown [3] that the description k = aT + 6T 3 holds very well in superconducting 
YBa 2 Cu 3 0j / (YBCO) samples below ~110 mK. 

Recently, there is a controversy in the high-T c community as to how best the k(T) data at 
low temperature are analyzed to separate K e from k p . In a paper reporting the low-temperature 
k(T) behavior in YBCO and La2_ x Sr x Cu04 (LSCO), Sutherland et al. argued [5] that in 
single-crystal samples of those materials, the specular reflection of phonons becomes important 
at temperatures below ~0.5 K, and hence the analysis of the low-temperature data should take 
this effect into account. The new analysis proposed in Ref. [5] was adopted in subsequent papers 
from the same group, such as Refs. [6j [91 [TO] , and has been playing crucial roles in deducing 
the conclusions of those papers. In this article, we point out that the key assumption behind 
the new analysis, that the phonons are specularly reflected at T up to ~0.5 K, is erroneous, by 
showing that the phonon mean free path is much shorter than the mean sample width in most 
of the proposed temperature range, which means that phonons are not even in the boundary 
scattering regime that is a prerequisite for the dominance of specular reflections. 

2. Specular Phonon Reflection? 

The standard way of analyzing the low-T thermal conductivity data k(T) in terms of k/T = 
a + bT 2 , which was first employed in the present context by Taillefer et al. [llj for both 
insulating and superconducting YBCO, has proved to be reliable for all cuprate systems studied 
so far [T2J [TOJ O EH H7] . The physical ground of this analysis is that the T dependence of k 
does not change further at lower T when the electrons are elastically scattered and the phonons 
are scattered by the boundary, and hence one can safely assume the formula k/T = a + bT 2 
to be valid down to T = K, which gives confidence in the T — > extrapolation; note that 
in the boundary scattering regime the phonons are usually diffusively scattered at the sample 
surface |18j and the phononic contribution k p is simply proportional to the specific heat, yielding 
k p = 6T 3 . However, Sutherland et al. had advocated in their paper [5] that the specular 
reflection of phonons causes a T dependence of k p weaker than T 3 in the boundary scattering 
regime, and that a fitting of the data for a wider temperature range to k/T = a + bT" -1 with 
a < 3 is better than the standard one. In the following, since the analysis of the n(T) data of 
underdoped YBCO has been the most controversial (HO [9] after the publication of Ref. [5], we 
focus our discussion on the case of YBCO. 

When the phonons are specularly reflected at the surface, k p becomes dependent on the 
averaged phonon wave length that changes with T, and empirically k p is expressed as bT a with 
a < 3 |18l I19j . Past studies on insulating materials unanimously found |18[ \19\ [20] that the 
phonon mean free path t becomes much longer than the mean width W of the sample when 
the specular reflection becomes important, as is naturally expected (W is usually taken to be 
2/y/n times the geometrical mean width w, see Ref. [111). Therefore, calculating I in YBCO 
and comparing it with W is a good way of judging whether the phonons are indeed specularly 
reflected. This can be rather easily done, because k p is equal to ^Cv£ where C = (3T S is 
the phonon specific heat and v is the averaged sound velocity, and both (3 and v are known 
for YBCO [11] . Figure 1 shows the temperature dependence of £ calculated with the data 
reported by Taillefer et al. [11] for y = (where W was specified) and with our own data [1] 
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Figure 1. Temperature dependence of the calculated phonon mean free path £ divided by the 
mean sample width W for YBCO. The data for y = 6.0 and 6.45 are taken from Refs. [TT] and 
[I], respectively. 

for y = 6.45 [21J. One can see that £ becomes saturating and comparable to W only at the 
lowest temperature, and it stays smaller than W in the temperature range where Sutherland et 
al. argued [HI E] that the specular reflection is taking place (0.07 - 0.5 K). Therefore, Fig. 1 
demonstrates that the assertion of Sutherland et al. is dubious and the specular reflection is not 
important in YBCO. A natural conclusion is that in YBCO at temperatures higher than ~120 
mK phonons are not in the boundary scattering regime |1 1|. [T5j [IB] and the phenomenological 
T a (a < 3) dependence comes from some additional scattering (such as scatterings off point 
defects, dislocations, twin boundaries, etc.) that only reduces £ [22]; such additional scattering 
would eventually die away at lower T as the boundary scattering regime is achieved, so the T a 
fitting should not be used for the T — > extrapolation. 

3. Validity of Extrapolations 

As a matter of fact, one must always bear in mind that, if an extrapolation is to be employed 
in the analysis, the formula used for the extrapolation should have a physical reasoning which 
assures that the same functional form holds throughout the range of the extrapolation. Otherwise, 
it is always possible that the true temperature dependence of the system changes below the 
experimental range (due, for example, to a change in the dominant scattering mechanism) and 
the extrapolation gives an inadequate estimate for the T = value. In the present case, when 
the assumption of the specular phonon reflection fails, there is no physical ground to expect the 
phenomenological T a dependence of k p to continue to T = 0, and hence the conclusion about 
the electronic contribution at T = loses its footing. 

4. Conclusion 

We have shown that in high-quality YBCO single crystals, the phonons are not likely to be 
specularly reflected at ~0.5 K because the estimated phonon mean free path is much shorter 
than the mean sample width at such relatively high temperature. Hence, the data analysis 
assuming the specular phonon reflection lacks its physical ground. When analyzing the low- 
temperature thermal conductivity data, one should avoid falling into pitfalls of fitting the data 
using some accidental power law and extrapolating it to zero temperature, when there is no 
particular reason to expect to its validity. 
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